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ABSTRACT 

Spectroscopic observations of the fast X-ray transient and superluminal jet source SAX J1819. 3-2525 (V4641 
Sgr) reveal a best fitting period of Pspect = 2. 8 1678 ±0.00056 days and a semiamplitude of = 21 1.0±3.1 km s"' . 
The optical mass function is f(M) = 2.74±O.12M0. We find a photometric period of Pphoto = 2.8 1730 ±0.00001 
days using a light curve measured from photographic plates. The folded light curve resembles an ellipsoidal 
light curve with two maxima of roughly equal height and two minima of unequal depth per orbital cycle. The 
secondary star is a late B-type star which has evolved off the main sequence. Using a moderate resolution spectrum 
{R = 7000) we measure T^ff = 1 0500 ± 200 K, log g = 3 .5 ± 0. 1 , and Vrot sin /= 1 23 ± 4 km s"' ( 1 cr errors). Assuming 
synchronous rotation, our measured value of the projected rotational velocity implies a mass ratio of 2 = Mi /M2 = 
1 .50 ± 0.08 (Icr). The lack of X-ray eclipses implies an upper limit to the inclination of / < 70.7°. On the other 
hand, the large amplitude of the folded light curve (« 0.5 mag) implies a large inclination (/ > 60°). Using the 
above mass function, mass ratio, and inclination range, the mass of the compact object is in the range 8.73 < 
Ml < II.7OM0 and the mass of the secondary star is in the range 5.49 < M2 < 8.14Mq (90% confidence). The 
mass of the compact object is well above the maximum mass of a stable neutron star and we conclude that V4641 
Sgr contains a black hole. The B-star secondary is by far the most massive, the hottest, and the most luminous 
secondary of the dynamically confirmed black hole X-ray transients. We find that the a-process elements nitrogen, 
oxygen, calcium, magnesium, and titanium may be overabundant in the secondary star by factors of two to 10 
times with respect to the sun. Finally, assuming E{B-V) = 0.32±0.10, we find a distance 7.40 < d < 12.31 
kpc (90% confidence). This large distance and the high proper motions observed for the radio counterpart make 
V4641 Sgr possibly the most superluminal galactic source known with an apparent expansion velocity of > 9.5c 
and a bulk Lorentz factor of F > 9.5, assuming the jets were ejected during during one of the bright X-ray flares 
observed with the Rossi X-ray Timing Explorer. 

Subject headings: binaries: spectroscopic — black hole physics — stars: individual (V4641 Sgr) — X-rays: stars 



Based on observations collected at the European Southern Observatory, Chile (program 65.H-0360) and the William Herschel Telescope operated on the island of 
La Palma by the Isaac Newton Group in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias 

^ Also at Astronomical Institute, Utrecht University 
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1. INTRODUCTION 

SAX J1819. 3-2525 was discovered as a relatively faint X- 
ray source independently with the Wide Field Cameras on Bep- 
poSAX on 1999 February 20 (in 't Zand et al. 1999, 2000) and 
with the Proportional Counter Array on the Rossi X-ray Tim- 
ing Explorer on 1999 February 18 (Markwardt, Swank, & Mar- 
shall 1999). The source had two rapid and bright X-ray flares 
around 1999 September 15. The first flare had a peak intensity 
of about 4.5 Crab in the 2-12 keV X-ray band (Smith, Levine, 
& Morgan 1999) and 5 Crab in 20-100 keV X-ray band (Mc- 
Collough & Finger 1999). The second peak, which came about 
0.8 days later, reached intensities of 12.2 and 8 Crab in the two 
bands, respectively. A few days after the giant flares the source 
was no longer detected in X-rays [e.g. Wijnands & van der Klis 
(2000)]. 

SAX J1819. 3-2525 was detected as a bright radio source 
(0.4 Jy) on 1999 September 16 (Hjellming et al. 2000), about 
16 hours after the second X-ray flare. The radio source was 
marginally resolved in this first observation, indicating high ex- 
pansion velocities. Thereafter, the source was detected at ra- 
dio wavelengths for about three weeks. Hjellming et al. (2000) 
modelled the radio images and the detailed radio light curves as 
a combination of the ejection of a relativistic, freely expanding 
jet and a subsequent ejection of a more slowly decaying, opti- 
cally thin jet segment. The radio observations did not resolve 
the moving components of the jets, so as a result one must make 
an assumption about when the jets were ejected in order infer 
the proper motion. If the jets were ejected during the brief 4.5 
Crab X-ray flare, the apparent proper motion is 0.22 arcseconds 
per day. Hjellming et al. (2000) concluded the most likely time 
of the jet ejection was during the rise of the 12 Crab X-ray flare, 
in which case the apparent proper motion would be 0.36 arcsec- 
onds per day. These high proper motions obviously imply large 
velocities; an apparent proper motion of > 0.22 arcseconds per 
day corresponds to an apparent velocity of > 1.28cd, where c 
is the velocity of light and d is the distance in kpc. 

SAX J1819. 3-2525 was already known as an optical variable 
star (Goranskij 1978), which was discovered 21 years before 
the first known X-ray activity. The source had a « 2 mag flare 
in 1978 June, and otherwise showed variabihty with an ampli- 
tude of up to « 1 mag about a mean level of B » 14.2 (Goranskij 
1990). In spite of differences in the classifications and variabil- 
ity behavior, Goranskij's variable was confused with GM Sgr, 
a long-period variable star discovered by Luyten (1927). As 
a result, many of the papers discussing the optical counterpart 
of SAX J1819.3-2525 refer to GM Sgr. Hazen et al. (1999) lo- 
cated the plate material used by Luyten and found that Luyten's 
variable GM Sgr is actually one arcminute south of Goranskij's 
variable star, which is the star at the precise radio coordinates of 
the X-ray transient. Goranskij 's variable has been given the new 
variable star designation V4641 Sgr, and we shall use this name 
for the remainder of the paper Goranskij (1990) assembled a 
quiescent light curve using 345 plates taken over a span of about 
30 years with the Crimean 40 cm astrograph. He noted a peri- 
odicity at a frequency of 2.7151 cycles per day, corresponding 
to a double-wave period of 0.7365483 days. However, owing to 
severe aliasing problems, Goranskij (1990) could not rule out 1 
day abases of his preferred period. We demonstrate below that 
the orbital period is 2.8713 days, which is an aUas of the period 
noted above (i.e. 1/(0.5 * 2.8173) 0.71 cycles per day). 

A network of optical observers had noted an increase in the 
variability of V4641 Sgr just prior to the giant X-ray flares 



(Kato et al. 1999). About six days before these flares the source 
showed a « 1 mag modulation with a period of « 2.5 days. The 
optical brightness peaked at V = 8.8 on 1999 September 15.8, 
and subsequently decayed rapidly to its mean quiescent level 
within two days. An optical spectrum obtained 1999 Septem- 
ber 16.25 showed very strong Balmer emission fines (Djorgov- 
ski et al. 1999). The Balmer lines made an abrupt transition 
from emission to absorption between 1999 September 17 and 
19 (Garcia, McCfintock, & Callanan 1999). A spechnm ob- 
tained 1999 September 30 (Wagner 1999) showed only strong 
Balmer absorption fines and interstellar features. At the time, 
this was quite confusing since the quiescent counterpart was 
thought to be an early K-type star based on the spectrum of 
GM Sgr shown in Downes et al. (1995). However, in hind- 
sight, the Balmer absorption line spectrum makes sense since 
we now know that Downes et al. (1995) did not observe V4641 
Sgr, which in fact contains a late B star. It now appears that the 
quiescent optical state was reached by 1999 September 19 (i.e. 
about four days after the X-rays were undetectable). 

It was initially thought that V4641 Sgr was relatively nearby 
based on two arguments. First, if the intrinsic velocity of the 
jet(s) observed in V4641 Sgr is similar to the velocities ob- 
served in GRO J1655-40 and GRS 1915H-105 (Mirabel & Ro- 
driguez 1999), then the large observed proper motion of the jet 
observed in V4641 Sgr would imply d w 500 pc. Second, a K- 
star in a « 0.7 day orbit with an apparent V magnitude of « 13.7 
would have a distance of « 1 kpc. We demonstrate below that 
d > 7.4 kpc at the 90% confidence level. Thus the apparent ve- 
locity of the 0"25 jet(s) seen several hours after the beginning 
of the X-ray event was > 9.5c, assuming the jets were ejected 
at the start of the giant X-ray flare. 

The galactic "microquasars" are excellent laboratories for the 
study of relativistic jets since they evolve orders of magnitude 
more quickly than the jets in quasars evolve (Mirabel & Ro- 
driguez 1999). V4641 Sgr is of special interest because of its 
very rapid X-ray flaring behavior and its extreme superluminal 
jets, and because the source is optically bright in quiescence 
(y « 13.7). In this paper we report the results of our spectro- 
scopic observations of V4641 Sgr. The observations and basic 
data reductions are summarized in Section 2. In Section 3 we 
establish the basic orbital parameters of the system. The mod- 
elling of a high resolution spectrum and the derivation of the 
secondary star properties are discussed in Section 4. We discuss 
the astrophysical parameters of the binary system in Section 5 
and discuss the impfications of our results in Section 6. The 
paper ends with a short summary in Section 7. 

2. OBSERVATIONS 

We obtained a total of 45 spectra of V4641 Sgr between 1999 
September 17 and 1999 October 16 using the FAST spectro- 
graph (Fabricant et al. 1998) on the 1.5 m telescope at the Fred 
L. Whipple Observatory (FLWO) on Mount Hopkins, Arizona. 
The 300 g rmn~' grating was used with a Loral 512 x 2688 
CCD. About half of the observations were taken with a 1"5 
wide slit, yielding a spectral resolution of ss 4.1 A FWHM. 
The other observations were taken with a 5."0 slit, yielding a 
spectral resolution of k. 6.0 A FWHM. The spectra cover the 
wavelength range 3650-6800 A. The exposure times were from 
three to ten minutes and the spectra were generally obtained un- 
der mostly clear skies with seeing between 1 and 2 arcseconds. 
The spectrograph sUt was rotated to maintain an approximate 
afigrmient with the parallactic angle, and a wavelength cafibra- 
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tion lamp was observed after each exposure. 

As we noted above, the source was active until about 1999 
September 18. All of the spectra from September 17 and 18 
show a broad Ha emission line, with a full width at zero in- 
tensity of up to « 5000 km s"' . In contrast, the higher Balmer 
lines were always in absorption. By September 19, the Ha line 
was also in absorption, and apparently the source was near or 
at its quiescent level. In the present paper we will discuss 14 
spectra obtained in quiescence (i.e. between 19 September and 
16 October) that were taken through the 1 ."5 wide slit. 

A total of forty-seven additional spectra of the source were 
obtained 2000 June 4-7 using the FORSl instrument on Antu, 
which is the first 8.2 m telescope at the European Southern Ob- 
servatory, Paranal. Thirty-seven of the spectra were taken with 
the 600B grism and a 0"7 wide slit; this combination gives 
a spectral resolution of 4 A FWHM and a wavelength cover- 
age of 3366-5751 A. The other ten spectra were taken with the 
600R grism and a 0"7 wide slit, yielding a spectral resolution 
of 3.2 A FWHM and a wavelength coverage of 5142-7283 A. 
The seeing was relatively poor for the VLT (« 1 ."5), and clouds 
were present on the night of June 5-6. The exposure times were 
from 1.5 to 4 minutes depending on the conditions. An atmo- 
spheric dispersion corrector was used, so the slit was kept at 
the default north-south direction. Following the standard pro- 
cedure at Paranal, the flat-field and wavelength calibration ex- 
posures were obtained during the daytime hours with the tele- 
scope pointed at the zenith. 

Ten additional spectra of the source were obtained 2000 June 
16 and 17 with the 3.5 m New Technology Telescope (NTT) at 
the European Southern Observatory, La Silla, using the blue 
arm of EMMl. The instrumental configuration consisted of 
the Tek 1024 x 1024 CCD, a I'.'O wide slit, and the grating 
#11, which is a high efficiency holographic grating with 3000 
"grooves" per mm (Wilson et al. 1991). The spectral resolu- 
tion is 0.43 A FWHM and the wavelength coverage is 3870.6- 
4032.8 A. The sUt was aligned with the parallactic angle, and 
a thorium-argon lamp was observed for wavelength calibration. 
The weather was clear June 16 and a total of eight 30-minute 
exposures were obtained. The conditions on June 17 were quite 
poor, and only two 30-minute exposures were taken. 

Finally, seven additional spectra were obtained 2000 August 
12 and 13 with the 4.2 m William Herschel Telescope (WHT) 
at the Observatorio de Roque de los Muchachos, la Palma, 
Spain. The instrumental configuration was the blue arm of the 
ISIS double beam spectrograph, the R1200B grating, and the 
EEV12 4096 x 2048 CCD. We used a 0."8 wide sUt, yielding a 
spectral resolution of 0.58 A FWHM and a coverage of 3680- 
4500 A. The slit was aligned with the parallactic angle, and 
a copper-argon lamp was observed for wavelength calibration. 
The weather was photometric, and the seeing was generaUy 1 ." 1 
or less. 

We reduced all of the spectra using IRAF^. The standard 
tasks were used to apply the bias and flat-field corrections, 
and to extract and wavelength-calibrate the spectra. The bright 
night sky emission line at » 5578 A was used to make small 
adjustments to the wavelength scales of the FLWO and VLT 
spectra. The shifts required to align this feature to a connmon 
wavelength (5578.0 A) were generaUy less than w 15 km s~^ 

^ IRAF is distributed by the National Optical Astronomy Observatories, which are 
cooperative agreement with the National Science Foundation. 
http://www.kusastro.kyoto-u.ac.jp/vsnet/Mail/obs26000/msg00925.html 



although the VLT spectra from the end of the night of June 7 
required shifts on the order of 80 km s~' . 

3. ORBITAL PERIOD AND SPECTROSCOPIC ELEMENTS 

3.1. Spectroscopic Period 

We measured the radial velocities of the secondary star us- 
ing the fxcor task within IRAF, which is an implementation of 
the technique of Tonry & Davis (1979). Synthetic spectra with 
Teff = 10500, log^ = 3.5, and Vrotsin; = 123 km s"' (see below) 
were used as templates for the five different sets of spectra (i.e. 
FLWO, VLT "blue", VLT "red", NTT, and WHT). The cross 
correlation regions covered the available wavelengths between 
3700 and 6847 A, excluding the interstellar Ca II H and K and 
Na D lines, the diffuse interstellar bands near 4330 and 5780 A, 
and a teUuric feature near 6280 A. The final velocities for the 
FLWO, the VLT blue, and WHT spectra were insensitive to the 
exact starting and ending points of the cross correlation regions 
and to the function used to find the centroid of the cross corre- 
lation peak. On the other hand, the velocities of the VLT red 
spectra and the NTT spectra were sensitive to these parameters 
since these spectra contain very few lines. We can estimate ap- 
proximate velocity corrections for the VLT red and NTT spec- 
tra by using observations of an AO template star observed both 
with the VLT and NTT. However, we find that including the cor- 
rected velocities did not improve the resulting spectroscopic pa- 
rameters since the corrected velocities all have relatively large 
uncertainties. We therefore excluded the VLT red and NTT ve- 
locities from the analysis presented below. 

To search for the spectroscopic period we computed a three- 
parameter sinusoid fit to the velocities for a range of trial pe- 
riods between and 5 days and recorded the values of the re- 
duced for the fits. The free parameters at each trial period 
are the velocity semiamplitude K2, the epoch of maximum ve- 
locity ro(spect), and the systemic velocity 7. The error bars on 
the individual velocities were scaled by a factor of 2.65 to yield 
= 1.0 at the minimum. Fig. la shows the resulting vs. P 
curve. The best fit is at a period ofP = 2.81678 ± 0.00056 days 
(1(7 error). The next best fit occurs at P = 2.699 days, where 
xi = 2.27. This alias period and the others are clearly ruled out 
by inspection of the folded velocity curves and we adopt the 
following spectroscopic elements: fspect = 

2.81678 ±0.00056 

days, /sTz = 211.0±3.1 km s-\ 7= 107.4±2.9 km &-\ and 
ro(spect) = HJD 2,451,442.523 ±0.052. The resulting optical 
mass function is then /(M) = 2.74 ± 0. 12M0 (Icr errors). Fig. 2 
shows the velocities and the best fitting sinusoid, and the spec- 
troscopic elements are listed in Table 1. 

3.2. Photometric Period 

Goranskij (1990) published a light curve of V4641 Sgr ob- 
tained from 345 plates taken with the Crimean 0.40 m astro- 
graph between about 1960 and 1990. After it was realized 
that V4641 Sgr was the optical counterpart of SAX J1819.2- 
2525, Goranskij reanalyzed the photographic data and made 
them publically available via the VSNet"*. We searched these 
data for periodicities using both the technique of Stellingwerf 
(1978) and the Lomb-Scargle Periodogram (Lomb 1976; Scar- 
gle 1982). We excluded the uncertain measurements and also 
data from 1978 June, when the source had an optical flare. Fig. 

operated by the Association of Universities for Research in Astronomy, Inc., under 
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Fig. 1. — (a): The .solid line .shows the reduced xi for three-parameter sinusoid fit to the 71 radial velocities as a function of the trial period in days. The 
minimum xi is at P = 2.8177 it 0.0005 days, (b): The solid line shows the Lomb-Scargle periodogram of the photometric light curve. Peaks above the dashed line 
are significant at the 3(t level. The peak with the most power is at P = 1.40865 days, which is half of the spectroscopic period, (c): The Stellingwerf (1978) 
statistic as a function of the trial period near 2.8 days. The minimum and hence the most coherent modulation is at P = 2.81730 ± 0.00001 days, (d): An expanded 
view of the curve in (a). The various aliases of the spectroscopic period are ruled out by the photometry. 



lb shows the Lomb-Scargle periodogram in the period range 
of to 5 days. Using a Monte Carlo procedure, we computed 
the power above which the peaks are significant at the 3a level. 
We made 10,000 simulated Ught curves having the same ob- 
servation times as the actual light curve and with magnitude 
points normally distributed about the mean magnitude. (By 
construction, the simulated light curves have the same variance 
as the actual light curve.) A Lomb-Scargle periodogram was 
computed for each simulated Ught curve and the power of the 
strongest peak was recorded. The 3<7 threshold was then de- 
termined from the cumulative distribution of the peak powers. 
This threshold is shown by the dashed line in Fig. lb. There 
are several peaks significant at the 3a level. The peak with the 
most power is at a trial period of P = 1 .40865 days (0.7099 cy- 
cles day"'), which is half of the spectroscopic period. All of 
the other significant peaks are aliases of the main peak (e.g. 
1/1-1/1.40865= 1/3.447, etc.) and their sidebands. Fig. Ic 
shows the 6 statistic of Stellingwerf (1978) plotted as a func- 
tion of the trial period in the range 2.6 < P < 3 . 1 days (we used 
20 phase bins of width 0.05 to compute the statistic — this 
combination yielded the smoothest minima). The minimum @ 
is at a period of P = 2.81730 days. Yearly aliases of this period 
are also evident. Fig. Id shows the vs. P curve from Fig. la 
plotted over the same period range as in Fig. Ic. The period of 
P = 2.81730 days is the only trial period where there are min- 
ima in both the statistic and in the xt ^ curve. Thus the 



photometric and spectroscopic periods of V4641 Sgr are nicely 
consistent with each other. 

We computed the uncertainty in the photometric period by 
using another Monte Carlo procedure. We adopted a mean 
photometric error of 0.1 mag (Goranskij, private communica- 
tion) and constructed 10,000 simulated data sets by assigning 
a new magnitude to each time using jnew,; =)'(+0.1(t, where 
C7 is a random Gaussian deviate with variance of unity and 
a zero mean. For each simulated data set we computed the 
Lomb-Scargle periodogram and the Stellingwerf O statistic and 
recorded the resulting derived periods. The uncertainty in the 
photometric period was taken to be the standard deviation of 
the distribution of best periods from these simulated data. We 
find P = 2.81730 ± 0.00001 days. Fig. 3 shows the fight curve 
folded on this period and binned into 20 phase bins. The light 
curve resembles an ellipsoidal light curve with two maxima of 
roughly equal heights and two minima of unequal depths per 
orbital cycle. The amplitude is quite large; the range is about 
0.5 mag from the deeper minimum to the maximum. 

Goranskij (1990) gives the time of minimum light 7()(photo) 
as HJD 2,447,707.454 (no error is quoted). We measured 
7o(photo) by folding the light curve on a dense grid of trial val- 
ues of ro(photo) and recording the of an ellipsoidal model 
fit (the solid line in Fig. 3). The result is ro(photo) = HJD 
2,447,707.4865 ±0.0038. The photometric phase of the time 
of maximum velocity ro(spect) is 0.75 ± 0.02, exactly as ex- 
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Fig. 2. — The radial velocities for V4641 Sgr folded on the spectroscopic period and phase and the best fitting sinusoid. See Table 1 for the parameters. Each 
point has been plotted twice for clarity. 




Phase 

Fig. 3. — The photometric light curve of V4641 Sgr folded using P = 2.81730 days and Tq = ro(photo)+P/2 and binned into 20 bins of equal size. Each point has 
been plotted twice for clarity. The error bars on the points represent the error of the mean. The solid curve is a i?-band ellipsoidal model with i = 70° ,2=1.6, and 
a grazing eclipse of the star by an accretion disk. The dash curve is the resulting hght curve when the disk is removed. 



pected for an ellipsoidal light curve. (Note that the phase zero 
point used in Fig. 3 is ro(photo)+0.5P to conform with our el- 
hpsoidal modelling code.) 

4. PARAMETERS FOR THE SECONDARY STAR 

V4641 Sgr is relatively bright, and as a result one can eas- 
ily obtain a spectrum with a relatively high signal-to-noise ra- 
tio (SNR) and high resolution. The eight NTT spectra from 
June 16 were Doppler shifted to zero velocity and combined 
to yield a "restframe" spectrum with a resolution of R = 9000 
and a SNR of about 50 per pixel. In the same way, we also 
made a restframe spectrum using the WHT data. The SNR of 
the WHT restframe spectrum is about 130 per pixel and is rea- 
sonably constant over much of the spectrum. We used synthetic 
spectra generated from solar metallicity Kurucz models^ to de- 

^ http://cfaku5.harvard.edu/ 



termine three basic parameters for the secondary star, namely 
its effective temperature Teff, surface gravity \ogg, and mean 
projected rotational velocity Vrotsini. The IDL program SYN- 
PLOT (Hubeny, Lanz, & Jeffrey 1994) was used to compute 
the detailed model spectra. S YNPLOT can generate an optical 
spectrum with arbitrary spectral resolution, wavelength sam- 
pling, and rotational broadening. [Currently SYNPLOT uses 
an analytic broadening kernel with a linear limb darkening co- 
efficient of e = 0.6 (Gray 1992)]. For the NTT spectrum the res- 
olution was taken to be 0.4333 A and the wavelength sampling 
was 0.1530 A, and for the WHT spectra those quantities were 
0.5802 A and 0.2231 A, respectively. We also have two syn- 
thetic spectra kindly provided by Peter Hauschildt computed 
using his PHOENIX code (Hauschildt & Baron 1999). These 
models are fuUy line-blanketed with magnesium and iron com- 



6 



\ 

\ 




110 115 120 125 130 135 

V^„tsin(i) (km s-') 

Fig. 4. — (a): The reduced as a function of the input value of Vrotsin; for the August 13 WHT restframe spectrum (solid Une) and August 12 WHT restframe 
spectrum (dashed Une), computed using the model spectra broadened by SYNPLOT. (b): The reduced as a function of the input value of Vmt sini for the August 
13 WHT restframe spectrum (solid line) and August 12 WHT restframe spectrum (dashed line), computed using the numerical broadening kernels. 



puted in non-LTE (NLTE). 

We used the technique outlined in Marsh, Robinson, & Wood 
(1994) to fit the model spectra to the observed ones. This fit- 
ting procedure has the advantage that it allows for the possi- 
bility of a (continuum) contribution from the accretion disk. 
The restframe spectra were normalized to their continuum fits. 
An input model spectrum characterized by the three parameters 
Teff, logg, and Vrotsin; is normalized to its continuum using the 
same procedure used for the observed spectra: it is scaled by 
a weight factor w, and subtracted from the observed spectrum. 
The scatter in the difference spectrum is measured by comput- 
ing the reduced from a low-order polynomial fit (after the in- 
terstellar lines were masked out of the fit). We varied T^tf, logg, 
and Vrotsini and looked for the "smoothest" difference spec- 
trum. The Kurucz grid has models for gravities of logg = 3.0, 
3.5, and 4.0 (cgs units) near the temperatures of interest (about 
10000 K). 

4.1. Effective Temperature and Spectral Classification 

It quickly became clear using the NTT spectrum that the 
gravity of the secondary star in V4641 Sgr is quite close to 
logg = 3.5 since the logg = 3.0 models have Balmer lines that 
are much too narrow and the logg = 4.0 models have Balmer 
Unes that are much too broad. We therefore fixed the gravity 
at logg = 3.5. The temperature is best constrained by the WHT 
spectrum. The He I line at 4024 A is somewhat sensitive to 
the temperature, and we find T^a = 10500 K. The statistical er- 
rors in the temperature and gravity are much smaller than the 



respective grid spacings (0.5 dex in logg and 250 K in Teff). 
No doubt there are systematic errors as well, but it is beyond 
our present ability to quantify these errors. We conservatively 
adopt (Tr = 200 K and aiogg = 0.1. The values of T^n and logg 
we find for V4641 Sgr are not too different from those of the 
B9 III MK spectral type standard 7 Lyr for which KiinzU et al. 
(1997) list Teff = 9970 ± 540 and logg = 3.50; the gravity was 
estimated from a measurement of the H/3 index (Balachandran 
et al. 1986). We fit synthetic spectra to the observed spectrum 
of 7 Lyr (obtained with the WHT) and find T^s = 10000 K, in 
good agreement with the previous values, and log g = 3 .0, some- 
what lower than the previous values. Based on this approximate 
similarity, we assign a B9111 classification to the secondary star 
in V4641 Sgr. We note that the spectral type of B91I1 derived 
here is somewhat earUer than the preUminary type of A2V we 
announced in lAU Circular 7440. The latter value was arrived 
at via a visual inspection of the Balmer lines in a low resolution 
spectrum. We believe the spectral type of B9111 derived here is 
much more reliable since we used an impartial fitting procedure 
on a spectrum with higher resolution and a higher SNR. 

4.2. Rotational Velocity of the Secondary 

Two other interesting facts also came to fight in this process. 
First, it was apparent that the rotational velocity of V464 1 Sgr is 
quite high (> 0.5A'2) since the atomic lines are quite broadened 
by rotation. Second, some of the metal lines (most noteably Mg 
II near 4781 A) are much stronger in V4641 Sgr than they are 
in the models. We therefore iterated the fits to determine the 
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Velocity (km s ') 

Fig. 5. — Top: The solid line is a numerical rotational broadening kernel which is the mean of 5 separate kernels matched to the August 13 WHT spectra, scaled 
by 125 km s"'. The dashed line is the analytic kernel used by SYNPLOT, scaled by 118 km s"'. Bottom: The solid line is a numerical rotational broadening 
kernel which is the mean of the two separate kernels matched to the August 12 WHT spectra, scaled by 130 km s~' . The dashed line is the analytic kernel used by 
SYNPLOT, scaled by 1 18 km s'^ . 



rotational velocity and the approximate abundances of Mg and 
other metals. 

We started first with the determination of the rotational veloc- 
ity. We computed model spectra for various values of Vrot sin i 
(using Jeff = 10500 K and log^ = 3.5) and computed the min- 
imum reduced of the polynomial fit to the difference spec- 
tra, masking out the few anomalous lines. It is well known 
that there is a systematic error in the measurement of Vrot sin; 
caused by the non-spherical shape of the Roche lobe and by 
the variations in the temperature and hmb darkening over the 
Roche surface. Therefore we made nightly restframe spectra 
using the two WHT spectra from August 12 and the five WHT 
spectra from August 13, thus insuring that the restframe spec- 
tra are composed of individual spectra with similarly distorted 
line profiles. Fig. 4a shows the reduced a function of 

the input Vrot sin / for the two restframe spectra. The minimum 
reduced x^ for the August 13 spectrum (which has the larger 
SNR) is for Vrot sin j = 118 km s~' . However, the curve is some- 
what noisy, and there is another dip near Vrot sin «' = 11 3 km s~^ . 
The corresponding curve for the August 12 restframe spectrum 
has its minimum at Vrot sin; =124 km s~', and another broad 
minimum near Vrot sin j = 1 30 km s~' . 

It is somewhat difficult to choose the correct value of Vrot sin i 
on the basis of these results. We repeated the exercise using 
phase-specific broadening kernels, computed using the ELC 
code (Orosz & Hauschildt 2000). The ELC code uses specific 
intensities from model atmosphere computations rather than 



from blackbodies and a one- or two-parameter Umb darkening 
law. Thus the distortions in the broadening kernel due to tidal 
distortions and to temperature and limb darkening variations 
are accounted for in detail. For the present problem we used 
the specific intensities from the Kurucz grid. We computed a 
broadening kernel for each of the two nights by averaging ker- 
nels made for each individual spectrum which went into the 
restframe average. We used a custom made IRAF routine to 
compute the convolution. Fig. 4b shows the resulting reduced 
X^ vs. Vrot sin J. The two curves are much "cleaner", and both 
have well-defined minima near Vrot sin; = 123 km s"'. Fig. 5 
shows the numerical kernels used. In the case of the August 13 
spectrum, the analytic kernel scaled by 118 km s"' has almost 
the same full width at half maximum as the slightly asymmet- 
ric numerical kernel, scaled by 123 km s"^ Thus, to first order, 
we would expect a correction from 1 18 to 123 km s"\ which is 
what we find. Similarly, the analytic kernel scaled by 130 km 
s"^ approximately matches the full width at half maximum of 
the numerical kernel for August 12, scaled by 123 km s~^ . Thus 
we adopt Vrot sin/ = 123 km s~'. The curves in Fig. 4b start to 
become steep at ^ 119 and 127 km s"\ so we adopt a Icr error 
of 4 km s"^ 

4.3. Abundance Anomalies 

After determining the rotational velocity we then examined 
in more detail the metal Unes that were unexpectedly strong. 
SYNPLOT can alter the abundance of any given element and 
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Fig. 6. — The results of our exploratory abundance analysis for V4641 Sgr. We show the solar abundance base Kurucz model, and the models with the abimdance 
patterns altered to approximately match the line profiles. The observed spectrum is from the WHT (resolution 0.58 A FWHM). The spectrum was rebinned by a 
factor of 2.6 for comparison with the Ca n line and the O I line. The interstaUar Ca n lines are labelled with "IS". 



compute an approximate line profile. The resulting line profile 
is reasonably accurate provided the abundance of the element 
in question is not altered by more than a factor of « 3 from 
its abundance in the input model atmosphere (I. Hubeny, pri- 
vate communication). Since we currently cannot compute a 
fuU grid of the fuUy line-blanketed Kurucz models with arbi- 
trary abundance patterns, our results should be taken as only 
indicative. We computed models with altered abundances and 
compared them to the observed V4641 Sgr spectrum. In search- 
ing for elements with possible abundance anomalies, we were 
guided by the work of Israelian et al (1999) who found that the 
abundances of the a-process elements nitrogen, oxygen, mag- 
nesium, silicon, and titanium in the F-star secondary of GRO 
J1655-40 are 6-10 times solar. We concentrated on these same 
elements in V4641 Sgr. As a check, we also computed mod- 
els for 7 Lyr with altered abundances. The results are given in 
Table 2 and in Figs. 6, 7, and 9. Details for each element are 
provided below. 

Taken at face value, the LTE Kurucz models indicate that 
the magnesium in V4641 Sgr is enhanced by about a factor 
of five to seven with respect to solar However, we cannot fit 
the magnesium lines in 7 Lyr using the LTE Kurucz models 
(Fig. 7), which is an indication that a NLTE treatment is needed. 
In Figure 8 we show the Kurucz models and the Hauschildt 
models for V4641 Sgr and 7 Lyr near the Mg II A4481.22 
line. In the case of V4641 Sgr (T^ff = 10500, log^ = 3.5), the 
NLTE treatment makes almost no difference to the strength of 
the A4481.22 line. In contrast, the NLTE treatment makes the 



A4481.22 fine much stronger in 7 Lyr (Jeff = 10000, logg = 3). 
We conclude the unusual strength of the A4481.22 Une in 
V4641 Sgr is due to an enhancement of magnesium of a fac- 
tor of five to seven times solar. 

The pair of lines at A4128.04 and 4130.58 due to Si II indi- 
cates roughly solar abundance in V4641 Sgr; but as with mag- 
nesium, these Si II lines are not well modelled in 7 Lyr. NLTE 
effects probably are important for Si II as well. 

Perhaps the most convincing case is for oxygen. The weak 
O I Une at A4368.2 observed in V4641 Sgr requires an oxy- 
gen enhancement of 3 times solar, whereas the solar abundance 
model for 7 Lyr fits the line profile in that star. For compari- 
son, Balachandran et al. (1986) measured an oxygen abundance 
in 7 Lyr of about 1.3 times solar. The unresolved O I triplet 
near A6454.4 appears in the VLT spectrum (resolution 3.2 A 
FWHM, Fig. 9). Its strength suggests the abundance is even a 
bit more than three times solar (we did not attempt any further 
iteration, given the low resolution of the VLT spectrum). 

The steUar Ca II fine at A3933.66 in V4641 Sgr was resolved 
from the interstellar line in the August 12 WHT spectrum. A 
calcium enhancement of twice solar is needed to fit the line. 
The Ca II line in 7 Lyr is fit by the solar abundance model. 

There are no strong titanium lines that are sensitive to the 
abundance in the WHT spectrum. The VLT spectrum (resolu- 
tion 4 A FWHM) contains two Ti II lines that appeared to be 
much stronger than the corresponding model lines (Fig. 9). A 
titanium enhancement of a factor of about 10 can explain the 
profiles of the A5226.5 and A5 129.2 Unes, but not the blend at 
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Fig. 7. — The abundance analysis for 7 Lyr. We show the solar abundance base Kurucz model (Tf-g = 10000 K, logg = 3.0), and Vrotsin; = 88 km s~', and the 
models with the abundance patterns found from the V4641 Sgr spectrum. The observed spectrum is from the WHT (resolution 0.58 A I^WHM). 
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Fig. 8.— The Mg n A4481.22 lines in V4641 Sgr (left) and 7 Lyr (right), observed with the WHT (resolution 0.58 A FWHM). The synthetic spectra made fi-om 
the LTE solar metaUicity Kurucz models are shown as the dashed lines and the HauschUdt model spectra (computed with magnesium in NLTE) are shown as the 
solid lines. The NLTE treatment of magnesium in 7 Lyr changes the strength of the line significantly, whereas there is little change in the line strength in V4641 
Sgr. 



« A5187. 

There are no strong nitrogen lines in the WHT spectrum. 
There is a line due to N I in the VLT spectrum (resolution 3.2 A 
FWHM) which is anomolous. A model with a nitrogen en- 
hancement of a factor of 10 solar gives a reasonably good fit 
to the profile. We also searched the WHT and VLT spectra for 



suitable sulphur lines, but we did not find any. 

To sunnmarize the results of our exploratory abundance anal- 
ysis, we find the following enhancements: nitrogen is w 10 
times solar (low resolution spectrum only); oxygen is « 3 times 
solar (low and high resolution spectra); magnesium is enhanced 
by about five to seven times solar; siUcon is roughly solar, but 
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Fig. 9. — (a): The restframe spectrum of V4641 Sgr observed with the VLT FORSl and grism 600B (spectral resolution 4 A I^WHM), the solar abundance model 
(Teff = 10500 K, logg = 3.5, and Vrotsini = 123 km s~'), and the Mg and Ti enhanced models, b: The restframe spectrum of V4641 Sgr observed with the VLT 
FORSl and grism 600R (spectral resolution 3.2 A FWHM), the solar abundance model (7;^ = 10500 K, logg = 3.5, and Vrotsini = 123 km s"'), and the N, O, and 
Mg enhanced models 



the LTE models may not be adequate; calcium is » 2 times so- 
lar (high resolution spectrum); and titanium is « 10 times solar 
(low resolution spectrum only). As we stated above, these re- 
sults should be treated with extreme caution. We do not claim 
that the Unes we examined are the most suitable ones for an 
abundance analysis, nor do we claim the solar metallicity Ku- 
rucz models are the most appropriate. Clearly we need better 
data (e.g. higher resolution spectra with lines from more than 
one ionization stage of the elements in question) and better 
models (e.g. NLTE treatment for all of the elements, etc.). 

Finally, we repeated the procedure to find the rotational ve- 
locity, using a model with altered oxygen, magnesium, and cal- 
cium abundances. We found no change in the value of Vrot sin i. 
Fig. 10 shows the normalized WHT spectrum, the best (normal- 
ized) model, and the difference between the two. The agree- 
ment between the model and the data is quite good. The accre- 
tion disk is faint; it contributes less than 1% of the light. For 
comparison, we also show the difference between the model for 
7 Lyr and the observed spectrum (spectrum e). 

5. ASTROPHYSICAL PARAMETERS 

Now that we have presented the observational facts concern- 
ing V4641 Sgr, we proceed to discuss the astrophysical impli- 
cations of these observations. 

5.1. Inclination Constraints 



A straightforward upper limit to the inclination comes from 
the lack of X-ray eclipses. As we noted earlier, the source was 
active as a relatively weak X-ray source between 1999 February 
and September and was detected on numerous occasions by the 
PC A instrument on RXTE. C. Markwardt (private communica- 
tion) informs us that the source was detected in X-rays several 
times near the inferior conjunction of the secondary star (i.e. a 
spectroscopic phase range of 0.75 ± 0.04). Assuming the X- 
rays come from a relatively small region R < 0.01 Rq centered 
on the compact object, one can show i < 70.7° for Q = 1 .5. 

In principle, one can determine the inclination by modelling 
the so-called "ellipsoidal" variations in the light curve of the 
B-star companion. The B-star fills its critical Roche lobe (there 
is ongoing mass transfer) and as such is quite distorted. As the 
distorted star moves about in its 2.8173 day orbit, the projected 
area on the sky changes, which in turn, gives rise to changes 
in the observed fiux. The net result is that one observes a light 
curve with two maxima per cycle (corresponding to the quadra- 
ture phases when the star is seen "side-on") and two minima 
per cycle (when the star is seen "end-on"). The amplitude of 
the light curve and the details of its shape is a function of the 
inclination. 

In practice, however, there are difficulties which can intro- 
duce systematic errors into the ellipsoidal modelling. For ex- 
ample, the accretion disk can contribute a significant amount 
of flux, thereby "diluting" the ellipsoidal Ught curve observed 
from the secondary. From an empirical point of view, the black 
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Fig. 10. — (a): The restframe WHT spectrum of V4641 Sgr, normalized to its continuum and offset upwards by 2.0 units, (b): A synthetic spectrum with 
T^ff = 10500 K, logg = 3.5, and Vrotsin/ = 123 km s~' offset by 1.5. (c): A synthetic spectrum with T^g = 10500 K, logg = 3.5, and Viotsin/ = 1 km s~^ not offset, 
(d): The observed spectrum minus the T^g = 10500 K, logg = 3.5, and Vrotsin/ = 123 km s"' model, no offset, (e). The observed 7 Lyr spectrum minus the best 
fitting model (J^s = 10000 K, logg = 3.0, and Vrot sini = 88 km s"'), offset downwards by 0.4 units. 



hole binaries with cool companions (G, K, or M-type compan- 
ions) are subject to the biggest systematic uncertainties since 
their light curves change from one observing run to the next 
(McCUntock & Remillard 1986; Orosz et al. 1996; Haswell 
et al. 1993; Webb et al. 2000). On the other hand, the two 
black hole binaries with hot companions (GRO J1655-40 with 
an F6III companion and 4U 1543-47 with an A2V compan- 
ion) have light curves which are quite stable, and it is apparent 
that the underlying ellipsoidal light curve of the secondary star 
dominates in these two cases. The light curves of GRO 11655- 
40 have been recently modelled and reasonably tight inchnation 
constraints (to within about 5 degrees at 90% confidence) have 
been obtained (Greene, Bailyn, & Orosz 2000). Since V4641 
Sgr contains a hot companion, it seems reasonable to expect 
that its light curve will be dominated by the elUpsoidal modu- 
lations from the secondary. 

The only well-sampled light curve we have access to at the 
moment is the photometric light curve assembled by Goranskij 
(1990). Unfortunately, the individual measurements are rela- 
tively imprecise (the mean estimated error is 0.1 mag), and the 
folded and binned light curve (Fig. 3) has relatively large error 
bars. Nevertheless, the ellipsoidal modulation is apparent in the 
binned light curve. For the moment we will limit ourselves to 
some representative models and put off the more detailed mod- 
elling until precise CCD Ught curves in several filter bandpasses 



become available. We used the ELC code (Orosz & Hauschildt 
2000) to compute the models using specific intensities from the 
Kurucz grid. The dashed line in Fig. 3 is a Johnson B-band 
model for a lobe-filling star (Tmean = 10500 K) with an inclina- 
tion of / = 70° and a mass ratio of Q= 1.5. The amplitude of 
the model is too small (x^ = 35. 196 for the 20 data points). The 
solid line in Fig. 3 is the same model, but with a faint accretion 
disk added (the outer radius is 98% of the Roche lobe radius). 
The depth of the minimum at phase 0.5 is increased by the par- 
tial eclipse of the B-star by the disk. In this case, the amplitude 
of the model is much closer to what is observed, and the model 
fit is improved (x^ = 25.39 for the 20 data points). The model 
disk contributes about 0.25% of the light in the B-band, and this 
would be compatible with the observed limits (Sec. 4). To find 
an approximate lower limit to the inclination, we adjusted the 
inclination of the model until the x^ of the fit increased by 9. 
The mass ratio and the radius of the disk were not adjusted. We 
find /> 60°. 

We adopt an inchnation range that is uniform in the inter- 
val 60° < i < 70.7°. A partial eclipse of the secondary star is 
needed to reproduce the relative depths of the minima. These 
conclusions (apart from the upper limit on i which is a function 
of Q only) should be treated with caution, given that the light 
curve is relatively uncertain. 
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5.2. Binary Parameters 
For a Roche lobe filling star it is straightforward to show that 

^2 fl V 2 / 
where Q is the mass ratio, Rli {Q) is the sphere-equivalent ra- 
dius of the Roche lobe, and Vrot = 2-kRlJP. We used the ELC 
code to numerically compute Ri^ (Q) for V^i = 123 ± 5 km s~' 
and = 211.0±3.1 km s"'; we find Q = 1.50±0.13 (90% 
confidence). 

We used a simple Monte Carlo code to compute other in- 
teresting binary parameters (e.g. the component masses, sec- 
ondary star radius, etc.) and their uncertainties using four 
input quantities, namely the period {P = 2.81730 ± 0.00001 
days, la), the velocity semiamplitude of the secondary star 
{K2 = 21 1.0 ±3.1 km s"\ Icr), the projected rotational veloc- 
ity of the secondary star (Viot sin; = 123±4kms"', Icr), and the 
inclination (60.0° < ; < 70.7°, uniform distribution). The nom- 
inal masses of the black hole and secondary star are 9.6Mq and 
6.5Mq, respectively; the precise masses and other key results 
are summarized in Table 3. 

6. DISCUSSION 

The mass of the compact object in V4641 Sgr (9.6Mq) is 
well above the maximum mass of a stable neutron star, which is 
usually taken to be « IMq (Chitre & Harde 1976). We there- 
fore conclude that V4641 Sgr contains a black hole. V4641 
Sgr is also interesting in several other respects. The B-star sec- 
ondary is by far the most massive, the hottest, and the most 
luminous secondary of the ten dynamically confirmed transient 
black hole binaries. The total mass of the system, \6.2Mq, is 
also largest among the dynanoicaUy confirmed black hole bina- 
ries. The next most massive system is V404 Cyg whose total 
mass is about I2.7M3 (Shahbaz et al. 1994). The other eight 
confirmed transient black hole binaries have total masses which 
are probably below lOM©, although the uncertainties are large 
in many cases (Bailyn et al. 1998, and cited references). 

The secondary star in V4641 Sgr appears to be in an un- 
usual evolutionary state. To help illustrate this, we plot in Fig. 
1 lab the radius of the secondary and the luminosity of the sec- 
ondary as a function of its assumed mass. Since the density 
of a Roche lobe filling star depends only on the orbital period 
of the binary to a good approximation, these curves are essen- 
tiidly independent of the inclination or the mass ratio of the 
binary. We also show the radius and luminosity of single stars, 
taken from the "Geneva" stellar evolution models (Schaller et 
al. 1992). For our adopted mass range of 5.49 < M2 < 8.14M0 
(90% confidence), the secondary star is both smaller in size and 
significantly underluminous compared to single stars of simi- 
lar masses. An underluminous secondary star usually indicates 
that the star went through a phase of thermally unstable mass 
transfer A more detailed discussion of the evolutionary history 
of this star is beyond the scope of the current paper. 

Two remarkable features of V4641 Sgr, which are discussed 
below, are a consequence of its great distance. For an as- 
sumed secondary star mass, we know precisely the correspond- 
ing radius. The corresponding luminosity is also precisely 
known since the temperature is well known from the spec- 
trum. The distance then follows from the apparent V magni- 
tude and the extinction. Photometry of the source was obtained 
from CTIO 2000 July 23.1 (UT). We find B = 14.00 ±0.02 

* http://www.kusastro.kyoto-u.ac.jp/vsnet/Mail/vsnet/msg01879.html 



and V = 13.68 ±0.02. The intrinsic colors of this star are all 
close to zero, hence E{B -V) = 0.32. The star is of course 
variable, so for now we adopt an uncertainty on the B and 
V magnitudes of 0.1 mag. Fig. 11c shows the distance as 
a function of the assumed secondary star mass. Assuming 
Teff = 10500 ±200 K and 5.49 < M2 < 8.14Mq (90% confi- 
dence), we adopt 7.40 <d< 12.31 kpc (90% confidence). 

As we noted earher, V4641 Sgr was marginally resolved in 
the first radio map made after the 12 crab X-ray flare. The 
proper motion of the jet is uncertain since it is not known when 
the jet was ejected. Hjellming et al. (2000) give three possi- 
bihties for the jet proper motion, based on known X-ray/radio 
correlations in other sources: 

• /^jet = 0.22 arcseconds per day for an ejection at the start 
of the 4.5 crab X-ray flare, which implies Vapp > 9.47c 
forJ > 7.4 kpc; 

• A*jet = 0.36 arcseconds per day for an ejection at the start 
of the 12 crab X-ray flare (this is their preferred value), 
which implies Vapp > 15.4c for d > 7.4 kpc; 

• Mjet =1-1 arcseconds per day for an ejection at the end 
of the 12 crab flare, which imphes Vapp > 47.7c for 

d > 7.4 kpc. 

For comparison, all of the other Galactic sources with relativis- 
tic jets have apparent velocities of less than 2c (Mirabel & Ro- 
driguez 1999). Indeed, this source might be characterized as a 
galactic "microblazar", since the jet velocity is similar to what 
is typcially seen in blazars (Jorstad et al. 2001). It is worth- 
while to point out that there are other interpretations of the ra- 
dio observations. For example, Chaty, et al. (2001) note that 
the long-lived radio remnant seen in Fig. 5c of Hjellming et al. 
(2000) did not show any movement between 1999 September 
16 and 24. This raises the possibility that the radio emission 
was due to an interaction of ejecta with the surrounding inter- 
stellar medium at a distance of « 0.25J AU from the source, 
where the distance d is measured in parsecs. At a distance 
of 9.6 kpc, an angular size of 0.25 arcseconds corresponds to 
a physical size of 2400 AU. The corresponding hght crossing 
time is 13.9 days, so the jet would have had to been ejected 
on « September 2 at the latest in order to have interacted with 
matter at a distance of 2400 AU. Optical observers coordinated 
through the VSNet noted possible activity as early as August 
8^, so a jet ejection as early as September 2 cannot be ruled out. 

The apparent proper motion of the approaching jet depends 
on the intrinsic jet velocity /? = v/c, jet angle 6, and the distance 
d (Mirabel & Rodriguez 1999) 



(3sin6 c 
l-f3cosed' 



(2) 



Fig. 12 shows a contour plot of as a function of 9 and /3 
for assumed distances of d = 7.4 kpc (our lower limit at 90% 
confidence) and d = 9.6 kpc (our preferred value). A proper 
motion of fia > 0.22 arcseconds per day requires /3 > 0.9945 
and 9 < 12° for d = 7.4 kpc. The bulk Lorentz factor in this 
case is F = (1 - /3^)"'/^ > 9.5. The combination of d = 9.6 kpc 
and Ha = 0.36 arcseconds per day gives (3 > 0.999 (F > 22.4) 
and ^ < 6°. In any case, the bulk motions were extremely rel- 
ativistic, and the jet in V4641 Sgr was nearly aligned with the 
line of sight. On the other hand, the large inclination indicates 
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Fig. 1 1. — Top: The radius of the secondary stai" as a function of its assumed mass. The points are from the single star Geneva stellar evolution models. Middle: 
Similar to the top, except the luminosity of the secondary star is displayed. Bottom: The distance of the source (right y-axis scale) and the peak isotropic X-ray 
liuninosity (left y-axis scale) as a function of the assumed secondary star mass for three different values of the color excess. The photometry of the source in 
quiescence suggests E{B-V) = 0.32. 
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Fig. 12. — A contour plot of the apparent proper motion in arcseconds per day of the approaching jet in V4641 Sgr as a function of the jet angle and the intrinsic 
jet velocity in units of c. The dash countours are for a distance of 7.4 kpc and the solid contours are for a distance of 9.6 kpc. 



that the orbital plane is nearly perpendicular to the Une of sight. 
If the jet is aligned with the spin vector of the black hole, then it 
seems likely that the spin vector of the black hole is misaligned 
(by a wide margin) with the orbital angular momentum vector. 
The inclination of the transient jet observed during the 1994 
outburst of GRO J1655-40 was about 85° (Hjellming & Rupen 
1995), whereas the orbital inclination is close to 70° (Orosz 
& Bailyn 1997; van der Hooft et al. 1998; Greene, Bailyn, & 



Orosz 2000), indicating a misalignment of about 15°. Thus if 
the inclination of the jet can be taken as the incUnation of the 
black hole spin vector, then misaligned black hole spin vectors 
may be a common feature of these systems. 

The peak isotropic X-ray luminosity observed during the 
large X-ray flare in the 2-10 keV band is Lmax = 5 x 
10^^(^/0.4 kpc)2 erg s"' (Smith, Levine, & Morgan 1999). Us- 
ing our adopted distance we find logL^ax = 39.46 ± 0.23 (90% 
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confidence) for the 2-10 keV band. The luminosity in the 10- 
100 keV band is comparable to the luminosity in the 2-10 keV 
band, so the bolometric luminosity will be a factor of « 2 
higher. Thus, rather than the main outburst being substantially 
sub-Eddington, as initially thought (in 't Zand et al. 2000), the 
outburst was actually super-Eddington, since logLedd = 39.10 
foramassof 9.61M0. The bright soft X-ray transient 4U 1543- 
47 had a peak X-ray flux (1.5-3.8 keV) that was at least three 
times the Eddington limit (Orosz et al. 1998), so the super- 
Eddington peak X-ray flux observed for V4641 Sgr is by no 
means unique. 

Fender & Kuulkers (2000) have recently shown that the peak 
X-ray flux Px observed in the outbursts of the black hole tran- 
sients is positively correlated with the peak observed radio flux 
Pr. If all black hole X-ray transients produce a mildly relativis- 
tic jet (five sources have been resolved at radio wavelengths), 
then the existence of the PxIPr correlation implies that on aver- 
age the bulk Lorenz factors in the jets are relatively small (e.g. 
r < 5) since very large values of F would result in many sys- 
tems with radio emission beamed out of the line of sight and 
hence unobservable (Fender & Kuulkers 2000). These authors 
point out that there is nothing unusual about V4641 Sgr in terms 
of its peak X-ray and peak radio flux. The radio emission ob- 
served in V4641 Sgr was strongly beamed, and the fact that 
V4641 Sgr sits on the Px/Pr correlation suggests that the X- 
rays may also have been beamed and that the peak outburst lu- 
minosity was not necessarily super-Eddington. X-ray beaming 
may also be responsible for the super-Eddington outburst of 4U 
1543-47 and for the apparently large numbers (twelve or more) 
of the binary sources in M82 and Cen A with luminosities of 
up to 10"*' ergs s~' (Prestwich 2000). These sources may be 
~ IOM0 black holes with beamed X-ray emission rather than 
K, IOOM0 black holes emitting isotropically at the Eddington 
limit. 

Finally, we note that we can make an independent estimate of 
the distance if we assume that the large observed radial veloc- 
ity of the binary (Table 1) is due entirely to differential galactic 
rotation. For a source at a galactic longitude of £ = 6.77 and 
galactic latitude ofb = -4.79 we find d>l kpc using the rota- 
tion curve given in Fich, Blitz, & Stark (1989) and the standard 
lAU rotation constants of R[) = 8.5 kpc and Oo = 220 km s"'. 
This lower hmit is in good agreement with the lower limit we 
derived above. The space motion of the binary relative to its 
local standard of rest Vcon- can be computed from the distance 
we derived above and the observed 7 velocity. However, Vcon is 
unfortunately rather uncertain owing to the binary's close prox- 
imity to the galactic center. The distribution of Vcon values com- 
puted using the simple Monte Carlo code mentioned above is 
roughly uniform in the range -80 < Vcorr < 80 km s~' . 

Although our basic picture of V4641 Sgr is secure, there 
are many useful follow-up observations that can be done. The 
photometric light curve has relatively large errors; and Ught 
curves with greater statistical precision and for several band- 
passes should be obtained. Our tentative conclusion is that the 
secondary star is partially echpsed by the disk. It should be easy 
to confirm this with CCD Ught curves. The derived properties 
of the binary (e.g. component masses, etc.) are somewhat sen- 
sitive to the value of Vrot sin i. Thus another measurement of the 
rotational velocity of the B-star secondary should be obtained, 
preferably near the time of the inferior conjunction of the B- 
star, when the systematic bias introduced by the tidal distortion 
of the secondary star is minimized. FinaUy, since V4641 Sgr is 



so bright, one can easily obtain a high resolution, high signal- 
to-noise spectrum and do a more detailed abundance analysis. 
This was done for the black hole binary GRO J1655-40 by 
Israelian et al (1999) who reported a large overabundance of 
oxygen, magnesium, silicon, and sulphur, which they attributed 
to the capture of supernova ejecta from the progenitor of the 
present-day black hole. We find that the a-process elements 
nitrogen, oxygen, calcium, and titanium appear to be overabun- 
dant with respect to solar in V4641 Sgr. However, we do not 
find the same levels of enrichment in V4641 Sgr that were 
found for GRO J1655-40. The amounts of the various elements 
synthesized in a supernova explosion depend on the mass of the 
progenitor star and on the energy of the explosion (Nomoto et 
al. 2000). Thus the differences in the levels of enrichment of 
the a-process elements in GRO J1655-40 and V4641 Sgr may 
reflect the differences in the supemovae progenitors. 

7. SUMMARY 

We have presented the results of our spectroscopic campaign 
on the fast X-ray transient and superluminal jet source V4641 
Sgr. The spectroscopic period is fspect = 2.81678 ± 0.00056 
days and the radial velocity semiamplitude is^C2 = 211.0±3.1 
km s"^ The optical mass function is f{M) = 2.1 A ± O.12M0. 
For the secondary star we measure T^ff = 10500 ± 200 K, \ogg = 
3.5 ±0.1, and Vrotsini = 123 ±4 km s"' (la errors). 

The photometric period measured from an archival photo- 
graphic light curve is Pphoto = 2.81730 ± 0.00001 days. The 
light curve folded on the photometric phase resembles an el- 
lipsoidal light curve. Modelling of this light curve indicates a 
high inclination angle of i > 60° . The lack of X-ray eclipses im- 
plies an upper limit to the inclination of i < 70.7°. Using these 
inclination Umits and the above spectroscopic parameters, we 
find a compact object mass in the range 8.73 < Mi < 11 JOMq 
(90% confidence). This mass range is well above the maximum 
mass of a stable neutron star and we conclude that V4641 Sgr 
contains a black hole. 

The secondary is a late B-type star which has evolved off the 
main sequence. It is in a peculiar evolutionary state since both 
its radius and its luminosity are much smaller than the corre- 
sponding values for single stars with a similar mass. 

Finally, we find a distance in the range 7.40<J<12.31 kpc 
(90% confidence), which is at least a factor of w 15 larger than 
the initially assumed distance of « 500 pc. The peak X-ray 
luminosity was super-Eddington, and the apparent expansion 
velocity of the radio jet was > 9.5c. 
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Table 1 

Orbital Parameters for V4641 Sgr. 



Parameter Value 



Orbital period, spectroscopic (days) 2.81678 ± 0.00056 

Orbital period, photometric (days) 2.81730 ± 0.00001 

velocity (km s-i) 211.0 ±3.1 

7 velocity (km s"') 107.4 ± 2.9 

7b, spectroscopic^ (HJD 2,451,000+) 442.523 ±0.052 

To, photometric''(HJD 2,447,000+) 707.4865 ± 0.0038 

Mass function (M©) 2.74 ± 0. 12 

yrotsinj(kms-i) 123 ±4 



Note. — All quoted uncertainties are la. 

"The time of the maximum radial velocity of the secondary 
star. 

''The time of the deeper photometric minimum, correspond- 
ing to the superior conjunction of the secondary star. 



Table 2 

Results for abundance analysis. 



Line 


Spectral 

resolulion (A) 


V4641 Sgr 

Wx(A) 


solar abundance 

model Wx (A) 


adopted model 

Wx (A) 


abundance 

(x solar) 


NIA6483.7 


3.20 


0.115 ±0.007 


0.011 


0.078 


10 


OI A4368.2 


0.58 


0.099 ±0.002 


0.033 


0.056 


3 


OI A6453.6 
O I A6454.4 
O I A6456.0 


3.20 
3.20 
3.20 


0.226 ±0.005 


0.097 


0.115 


3 


Mgll A4481.2 


0.58 


0.455 ±0.010 


0.310 


0.476 


7 


Call A3933.7 


0.58 


0.573 ±0.010 


0.424 


0.550 


2 


Till A5 129.2 


4.00 


0.085 ±0.005 


0.019 


0.066 


10 


Till A5 185.9 
Till A5 188. 7 


4.00 
4.00 


0.083 ±0.012 


0.129 


0.247 


1 


Ti 11 A5226.5 


4.00 


0.1 10 ±0.005 


0.076 


0.128 


10 



Note. — 



All quoted uncertainties are la. 
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Table 3 

ASTROPHYSICAL PARAMETERS FOR V4641 SGR. 

Parameter Value 

Black hole mass (Mq) 9.61!^;^^ 

Secondary star mass (Mq) 6.53;|;} o3 

Total mass (Mq) 16. 19!^^^ 

Mass ratio 1.50 ±0.13 

Orbital separation (Rq) 21 .33!j;^f 

Secondary star radius (Rq) 7.47;I:q47 

Secondary star luminosity (Lq)^ 610;!;[^ 

Distance (kpc)'' 9.59!|]^ 

logZ,, (erg s-'r 39A6^il 



Note. — The quoted errors are all 90% confi- 
dence. 

^Assuming T^s = 10500 ± 200 K, la. 

''Assuming E(B-V) = 0.32 ±0.10, la and 
Av = 3.lE{B-V). 

'^The peak X-ray luminosity (2-10 keV), given 
by 5 X l0^\d/4 kpcf erg s"', (Smith, Levine, & 
Morgan 1999). 



